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Abstract 



, For a network of interconnected nonlinear dynamical systems an adaptive leader- 

Q I follower output feedback synchronization problem is considered. The proposed struc- 

ture of decentralized controller and adaptation algorithm is based on speed-gradient 
and passivity. Sufficient conditions of synchronization for nonidentical nodes are 
established. An example of synchronization of the network of nonidentical Chua 
systems is analyzed. The main contribution of the paper is adaptive controller de- 
sign and analysis under conditions of incomplete measurements, incomplete control 
^ , and uncertainty. 

00 

^ : 1 Introduction 

^ . Adaptive synchronization of netv\rorked dynamical systems has attracted a grov^ing inter- 

0\ ! est during recent years [1-4]. It is motivated by a broad area of potential applications: 
formation control, cooperative control, control of povv^er netv^^orks, communication net- 
v\rorks, production netvi^orks, etc. Existing v^orks [1-4] and others are dealing vi^ith full 
^ ■ state feedback and linear interconnections. The solutions are based on Lyapunov functions 
formed as sum of Lyapunov functions for local subsystems. As for adaptive control algo- 
rithms they are based on either local (decentralized [5-12]) or nearest neighbor (described 
by an information graph [13-16]) strategies. 

Despite a great interest in control of network, only a restricted class of them is currently 
solved. E.g. in existing papers mainly linear models of subsystems are considered [13,14]. 
In nonlinear case only passive or passifiable systems are studied and control is organized 
according to information graph, i.e. not completely decentralized [15,16]. Availability of 
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the whole state vector for measurement as well as appearance of control in all equations for 
all nodes is assumed in decentralized stability and synchronization problems [1-4]. Power- 
ful passivity based approaches are not developed for adaptive synchronization problems. 

In this paper we consider the problem of master-slave (leader-follower) synchronization 
in a network of nonidentical systems in Lurie form where system models can be split 
into linear and nonlinear parts. Case of identical nodes is studied in [17]. Linearity of 
interconnections is not assumed; links between subsystems can also be nonlinear. In the 
contrary to known works on adaptive synchronization of networks, see [3,4], only some 
output function is available and control appears only in a part of the system equations. 
It is also assumed that some plant parameters are unknown. The leader subsystem is 
assumed to be isolated and the control objective is to approach the trajectory of the leader 
subsystem by all other ones under conditions of uncertainty. Interconnection functions 
are assumed to be Lipschitz continuous. 

The results of [11,12] are employed to solve the posed problem. Adaptation algorithm 
is designed by the speed-gradient method. It is shown that the control goal is achieved 
under leader passivity condition, if the interconnection strengths satisfy some inequalities. 

The results are illustrated by example of synchronization in network of nonidentical 
Chua curcuits. 

2 Auxiliary results 

2.1 Yakubovich-Kalman Lemma 

We need Yakubovich-Kalman Lemma in following form, see [18]. 

Lemma 1 Let A,B,C be n x n,n x m,n x m real matrices and u G M™, xi.^) = 
C'^{sln — A)^^B ^i&VikB = m. Then the following statements are equivalent: 

1 ) there exists matrix H = > such that 

HA + A'^H<0,HB = C; (1) 

2) polinomial det(s/„ — A) is Hurwitz and following frequency domain conditions hold 

Reu^'x(iuj)u > 0, lim u'^ Reu^Y(iu)u > 

for allu e M\ M G W,u^ 0. 

2.2 Speed gradient algorithm in decentralized control 

In order to present the main syncronization result of this paper we need to formulate 
problem statement of decentralized control and Theorem 2.18 from [11] which can also be 
derived from Theorem 7.6 from [12]. 

Consideill] a system S consisting of d interconnected subsystems Si, dynamics of each 
being described by the following equation: 

Xi = Fi{xi,Ti,t) + hi{x,T,t), i = l,...,d, (2) 

^ In this paper norms are Euclidean, col(a;i, . . . ,Xd) stands for column vector with components con- 
sisting of components of Xi,i = 1, . . . , d. 



where xi G M"' - state vector, Tj G M™^ - vector of inputs (tunable parameters) of 
subsystem, x = col(xi, . . . , Xd) G M", r = col(ri, . . . , r^) G - aggregate state and input 
vectors of system n = J^^i, m = Y^f^i- Vector-function Fj(-) describes local dynamics 
of subsystem Si^ and vectors hi{-) describe interconnection between subsystems. 
Let Qi{xi,t) > 0,i = 1, . . . ,d he local goal functions and let the control goal be: 

Urn Qi(xi,t) = 0, i = l,...,d. (3) 

t — 5-00 

For all i = 1, . . . ,d we assume existence of smooth vector functions x*{t) such that 
Qi{x*{t),t) = 0, i.e. X* = aigmin^^ Qi{xi,t). Decentralized speed-gradient algorithm is 
introduced as follows: 

fi = -TiVr,uJi{xi,Ti,t), i = l,...,d, (4) 

where 

Ti = TJ > 0, rrii X rrii - matrix. 

Theorem 1 Suppose the following assumptions hold for the system S: 

1. Functions Fj(-) are continuous in Xj, t, continuously differentiable in and locally 
bounded in t > 0; functions Qi{xi,t) are uniformly continuous in second argument 
for all Xi in bounded set, functions uJi{xi, Ti, t) are convex in r^; there exist constant 
vectors r* G M™' and scalar monotonically increasing functions tii{Qi), Pi{Qi) such 
that Kj(0) = pi(0) = 0,limQ^_+oo/«i(<5i) = +00 

uji{xi, r*, t) < -pi{Qi{xi, t)), (5) 

and Qi{xi, t) > Ki(^\\Xi X 

m\\)- 

2. functions hi{x,T,t) are continuous and satisfy the following inequalities 

d 

\Vx,Qiixi,tyhi{x,T,t)\ <^Hijpj{Qj{xj,t)), (6) 

where matrix M — / is Hurwitz, M = {pij}, Pij > 0, / is identity matrix. 

Then system (IS]),^ is globally asymptotically stable in variables Xi — x*(t), all trajectories 
are bounded on i G [0, +00) and satisfy (l3|). 

3 Main result 

3.1 Problem statement. Adaptive controller structure 

Let the leader subsystem be described by the equation 

T = Alx + BLiu + Mv)). y = c^x, (7) 



where x G M" - state, y E - measurement, u{t) E M} is control that specified in 
advance, t/^q : — >^ - internal nonlinearity. Let Al, Bl, C and ^o(') be known and not 
depending on the vector of unknown parameters ^ G H, where S is known set. 

Consider a network S of d interconnected subsystems Si, i = l,...,d,d G M. Let 
subsystem Si be described by the following equation 

d 

Xi = AiXi + BiUi + BL^poiVi) + (^ij'^iji^i - ^j)^ 

j=i (8) 

Vi = C^Xi, i = 1, . . . , d, 

where Xi G M", Ui G IS}, aij G M}, iji G M'. Functions z = 1, . . . , (i, j = 1, . . . , d, de- 

scribe interconnections between subsystems. We assume ipa = (0,0,..., 0)^, i = 1, . . . ,d. 
Let matrices Ai, Bi and functions (Pij{-), i = I, . . . ,d,j = 1, . . . ,d, depend on the vector 
of unknown parameters ^ G H. 

Network model (IHl) can describe, for example, interconnected electrical generators [19]. 

Let the control goal be specified as convergence of all subsystems and the leader 
trajectories: 

\im {xi{t) -x{t)) = 0, i = l,...,d. (9) 

t — > + CXD 

The adaptive synchronization problem is to find a decentralized controller 
Ui = Ui{yi,u,t) ensuring the goal (I9|) for all values of unknown plant parameters. 

Denote ai{t) = col{yi(t) ,u{t)) . Let the main loop of the adaptive system be specified 
as set of linear tunable local control laws: 

Ui{t) =ri{tya,{t), z = l,...,ci, (10) 

where Ti{t) G R'"'"^,i = l,...,d are tunable parameters. By applying speed-gradient 
method [12] it is easy to derive the following adaptation law: 

T'i = -g^iVi - y)Tia^{t),i = 1, . . . , (11) 

where = > - (/ + 1) x (/ + 1) matrices, g G M'. 

3.2 Synchronization conditions 

Introduce the following definition. 

Definition 1 Let G G M}. Function /: — is called G-monotonically decreasing if 
inequality (x — y^G (/(x) — f{y)) < holds for all x,y eM.K 

Remark 1. Apparently, for / = 1,G = 1 G-monotonical decrease of the function / 
is equivalent to incremental passivity [21] of the static system with characteristics (— /). 
Definition [T] is easily extended to dynamical systems with the state vector x G M", input 
u G and output y G M'. It corresponds to existence of a smooth function y(xi,X2) 
satisfying an integral inequality 

V{xi, X2) < {u2 - uiyG{y2 - yi). 



The corresponding property can be called incremental G-passivity by analogy with [20]. 



Consider real matrices H = > 0, (7 of size nxn,lxl correspondingly and a number 
p > such that: 

HAl + AlH < -pH, HBl = Cg. (12) 

Denote A,, = \max{H) / \min{H) condition number of matrix iJ, where Xmax{H) , Xmin{H) 
are maximum and minimum eigenvalues of matrix H. 

For analysis of the system dynamics the following assumptions are made. 

Al) The functions (pij{-),i = 1, . . . ,d, j = 1, . . . ,d are globally Lipschitz: 

\\fij{x) - ipij{x')\\ < Lij\\x - x'W, Lij > 0. 

The function V'o(') is such that the unique existence of solutions of ((71) holds. 

A2) (Matching conditions, [22]) For each ^ G S there exist vectors z/j = z/j(^) G M} and 
numbers 9i = 9i{^) > such that for z = 1, . . . , ci 

Al = A + B,uJC\ Bl = e,B,. (13) 

Denote x{^) = C'^i.sln — Al)^^Bl. For the case when matrix Al is Hurwitz intro- 
duce notation p^, for stability degree of the function's g'^x{s) denominator, i.e. = 
v[im.k=i,...,n I R-eAfc(y4i^)| where Xk^Ai) are eigenvalues of A^. 

Theorem 2 Let B^ 7^ 0, matrix A^ he Hurwitz and for some g E the following 
frequency domain conditions hold: 

Re 5f'^x(^^^) > 0, Yim^u'^Reg'^xii^) > ^ (14) 

for all uj G Then there exist H = H'^ > 0,p > such that relations ( fT2ll hold. 

Let for all ^ G S Assumptions Al, A2 hold, function V^o(") be (^-monotonically decreas- 
ing, and following inequalities hold 

d 

^|aijLij|<7 z = l,...,(i, (15) 

i=i 

where 7 = p*/(4(iA*), A* is condition number of matrix H. 

Then for all ^ G S, z = 1, . . . , ci adaptive controller (fTOl) . (fTTl) ensures achievement of the 
goal (jni) and boundedness of functions Oiit) on [0, cxd) for all solutions of the closed-loop 
system ©, (i), 

Proof. Let's apply Lemma [TJ Note that in our case m = 1, i. e. u is scalar. Let's 
choose Cg instead of C in ([I]). Then statement of the Lemma[I]and conditions of Theorem 
[2] ensure existence of matrix H = H'^ > such that 

HAl + AlH < 0, HBl = Cg. 

Now we can conclude that there exists number p > such that the following is true: 

HAl + AlH < -pH, HBl = Cg. (16) 

Denoting introduce auxiliary error subsystems: 

d 

Zi =AiXi + BiUi + BLipoivi) + J2 oiij^ijixi -Xj)- {Alx + Bl{u + ^oiy))) , 

3=1 (17) 

Vi =C Zi., i = 1, . . . , d, 



here we choose Ui(t) same as in (fTOll . 

Let us choose following goal functions Qi{zi) = ^zfHzi, and apply Theorem [TJ We 
need to evaluate the derivative trajectories of Qi{zi) along trajectories of isolated (i.e. 
without interconnections) auxiliary subsystems f flTl) : 



Ui{xi,x,Ti) = zjH[AiXi + BiTj{t)ai{t) + BLi^oiVi)- Alx- BL{u + My))]- (18) 
Denote r* = col(z/j, 6'j), z = 1, . . . , d. By taking Tj = r*, i = 1, . . . , rf, we obtain 

uji{xi, X, T*) =zjH[AiXi + BiiviC^Xi + 6'iu) + Bi^poiyi) - Alx - Blu - Bi^Q^y)] = 
zJH[Alx, + Blu + Bl^iPq^v,) - ^Poiy)) - ALXi - Blu] = 

zjH[ALZ, + BLiMy^)~My))]- 

Further, for i = 1, . . . ,d 

zjHBLiMvi) - My)) = ^JCgiMvi) - My)) = (yi - WaiMyi) - My)) < o. 

The last inequality holds because V'o(") is (y'-monotonically decreasing. So 

Ui{xi,x, T*) < ^zJ{HAl + AlH)zi. 
Taking into account (fTGll we conclude 

Ui{Xi,X,T*) < -pQi{Zi). 

By taking Pi{Q) = p ■ Q we ensure that ([5]) holds for i = 1, . . . ,d. Other conditions from 
the first part of Theorem [U hold, since the right hand side of the system f flTl) and function 
Qi{zi) are continuous in Zi functions not depending in t for any i = 1, . . . ,d. Convexity 
condition is valid since the right hand side of (fTSl) is linear in r^. 
The interconnection condition (l6|) in our case reads: 

d d 

l^zM^i)^ J2 (^ijfiji^i - %)l < J2 f^ijp ■ Qi^j)^ (19) 

where i = 1, . . . ,d, and matrix M — I should be Hurwitz (M = {pij}, Pij > 0). 

For the case d = 1 we can take pu = 0.5 and last inequality will be satisfied. Let's 
consider case d > 1. 

For i = 1, . . . ,d rewrite (fT9l) as follows: 



z, 



H J2 o^iifij (zi- Zj)\<-Y, f^ij zjHzj . (20) 



Evaluate the left-hand side of (EHl 

d 

zjH J2 c^ij^ijizi - Zj) <] 



< J2 \zjHaij(fij{zi - Zj)\< 

d d 

j=i i=i 



\Zif + \\Zi\\ ■ \\Zj\\), 



ior i = 1, . . . ,d. Then for i = 1, . . . ,d evaluate lower bound of the right-hand side of (l20ll : 

It is seen that for i = 1, . . . , d to ensure ([6]) it is sufficient to impose an inequality 

d p d 

\aijLij\ ■ Xmax{H) ■ {\\Zif + \\Zi\\ ■ \\Zj\\) < - J2f^ij^min{H)\\Zjf, 
3=1 ^ i=i 



or 



d „ \ I TJ\ d 



■ {\\z^f+ Ikill ■ \\zj\\) < ^ ^"'"^.^1 IZ/^»ilkjir. i = l,...,d. (21) 



Denote ( = p/{amax ■ 2A*), where 

d 

i:l<i<d 



amax = max V \aijLij \ . 



Noting that p in ( fT2l) can be chosen arbitrarily close to and taking into account ( fTSl l 
we can conclude that 

C>2d. 

The left-hand side of (EB: 



d / d \ / d 

J2\aijLij\i\\zif + \\zi\\-\\zj\\) < ■ Ildkif + Ikill ■ IkjII) I < 

1 / ' \ 

-amax ISdWziW^ + Y2j\zj\\'^ ] , i = 

Thus, if following inequality holds then ((61) is ensured: 

d d 

3d\\z,f + Y.\\zjf<2C-J2f^ij\\^jf^ 2=l,...,c?. (22) 



Introduce matrix M = {fiij} as follows 



M 



( P'll P'12 ■■■ P'ld\ 
P-21 P-22 ■ ■ ■ P'2d 

XP'dl P'd2 ■ ■ ■ P'ddJ 



P'ij 



^(3rf+l), i=j; 



1 

2C' 



« 7^ J- 



Such choice of M ensures (I22]) . 

Note that M is symmetric. If matrix / — M is positive definite then M — / is Hurwitz. 
Diagonal elements of / — M are positive since d > 1 and C > 2c?. By taking into account 
that 

1 - ^(3d+l) - (d-1)— > 
2C 2C 



and applying Gershgorin circle Theorem we conclude that M — I is positive definite. 
Thus, statement of the Theorem [2] follows from Theorem [TJ □ 

Remark 2. The value of 7 can be evaluated by solving LMI (flGll by means of one of 
existing software package. 

Remark 3. By interconnections graph of network S we can consider directed graph 
which is a pair of two sets: a set of nodes and a set of arcs. Cardinality of a set of nodes 
is d] i-th node is associated with subsystem Si for any i = 1, . . . ,d. We say that arc from 
z-th node to j-th node belongs to the set of arcs if (Pij{-) is not zero function. By weighted 
in-degree of i-th node we define following number: J2'j=i WijLij \ ■ If each nonzero addend 
from last sum is equal to 1 then introduced definition of weighted in-degree of the node 
coincides with the definition of in-degree of digraph's node. Thus the inequality f fTSD can 
be interpreted as follows: weighted in-degree of each node of interconnections graph must 
be less than 7. 



4 Example. Network of Chua circuits 
4.1 System description and theoretical study 

Chua circuit is a well known example of simple nonlinear system possessing complex 
chaotic behavior [23]. Its trajectories are unstable and it is represented in the Lurie form. 
Let us apply our results to synchronization with leader subsystem in the network of five 
interconnected nonidentical Chua systems. 

Let mo = —8/7, mi = — 5/7,p = 15.6, g = 30, 6 = 1 and g = 1- 

Let the leader subsystem be described by the equation 

i = AlX + BL{u + i)Q{y)), y = C'^x, 

where x e is state vector of the system, y G is output available for measurement, u is 
scalar control variable, ■?/'o(y) = pv{y) /b, where v{x) = — 0.5(mo— mi)(|x+l| — |x — 1|— 2x). 
Further, let 

/-I 
Al= 1 -1 
V -q 

Bl = col(6, 0, 0), C = col(l, 0, 0). 

Transfer function = C^isI - AlY^Bl = (5^ + 5 + ?>Q)/{s^ + 2s'^ + 31s + 30). It 
is seen from the Nyquist plot of x(ii-i^), Vtu G M^, presented on Fig. [U that first frequency 
domain inequality of (fT4l ) holds. The second frequency domain inequality of ( fT4ll also 
holds since relative degree of is equal to one and highest coefficient of its numerator 
is positive. 

Obviously is (^-monotonically decreasing. 

Let subsystem S'j for i = 1, . . . , 5 be described by (JHl) with Uj, G M^. By choosing 
(z/i, z/2, t's, 1^4, t's) = (3, 1, 4, 1, 5), 6'j = = 1,...,5 and using (fT3l ) we obtain matri- 

ces Ai, Bi for i = 1,...,5, which are not equal, i.e. nodes are nonidentical. Denote 
^ij = ^ij{xi - Xj),i = l,...,5,j = 1,...,5. Let v^m, V^25, V^42, V^45, V^52, V^sa, be equal 



Figure 1: Nyquist plot of x{i^),^ ^ 



to (0, 0, 0)^. Further, let 



fl2 = 


(sm(xii -X2i),0,0)'^, 


V^is = 


{0,Xi2 - 0:32,0)'^, 


fl5 = 


(0,0,sin(xi3 - X53)y, 


f21 = 


(X2l — Xii, 0, X23 — 2:13)^ 


^23 = 


(0,Sm(x22 - X32),0)^, 


= 


(0,X22 - a;42,0)^, 


f31 = 


(sm(x3i -xii),0,0)^, 


V534 = 


(sin(x3i -X4i),0,0)^, 


^^35 = 


(Xsi — X51, X32 — X52, X33 — ^53)'^, 


V941 = 


(0,Sm(x42 - Xi2),0)^, 


<^43 = 


(sin(x4i - X3i),0,0)^, 


V^Sl = 


(X51 -a;ii, 0,0:53 -xisY 


V^54 = 


(0,X52 - X42,0)^. 







Lipschitz constants of all ipij are equal to 1. 

It follows from Theorem [2] that decentralized adaptive control (fTOll provides synchro- 
nization goal ([9]) if for alH = 1, . . . , 5 inequality Z]j=i \c(ij\ < 7 holds, i.e. if interconnec- 
tions are sufficiently weak. 



4.2 Simulation results 



Consider following control of leader subsystem u = j [(—(1 + ^«o)p + + PX2] ■ Such u 



ensures chaotic behavior of leader subsystem. Let us put Fj 
identity matrix, and 



, d, where I 



xi(0) = 0.5, X2(0) 
xi(0) = (7, 14,0.4)^ 
X3(0) = (1,-1,4.5)^ 
X5(0) = (2,8,15). 



0, X3(0) = 0, 
X2(0) = (0,4,4)- 
X4(0) = (3,-4,0.2) 



Denote by a 5 x 5 matrix with element ctj,- lying in the i-th row and the j'-th column. 



1, . . . , 5, and 



a 



( 0.0051 0.1395 0.1676\ 

0.0662 0.0921 0.0065 

0.2013 0.2271 0.1430 

0.0907 0.0675 

^0.0663 0.2773 / 



(A): (B): 




Figure 2: (A): Phase portrait of leader subsystem, (B): ||2;j||, (C): Ui = Ui — u,i = 
1,...,5. 

Let us choose adaptive control Ui,i = 1, . . . , 5 as in (fTOl) and apply Theorem [21 If we take 
a = a, then simulation shows that ||zj|| — > 0, z = 1, . . . , 5, i.e. synchronization is achieved: 
all state vectors of nonidentical nodes converge to the state vector of the leader subsystem, 
see Fig. [2]-(B). Phase portrait of the leader subsystem, \\zi\\,u = Ui — u,i = 1, . . . ,5 found 
by 40 sec. simulation are shown on Fig. [2l 

5 Conclusions 

In contrast to a large number of previous results, we obtained synchronization condi- 
tions for networks consisting of nonidentical nonlinear systems with incomplete mea- 
surement, incomplete control, incomplete information about system parameters and cou- 
pling. The design of the control algorithm providing synchronization property is based 
on speed- gradient method [12], while derivation of synchronizability conditions is based 
on Yakubovich-Kalman lemma and result presented in [11]. 
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